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Abstract Small amplitude oscillations are a commonly observed feature in prominences/fi- 
laments. These oscillations appear to be of local nature, are associated to the fine structure of 
prominence plasmas, and simultaneous flows and counterflows are also present. The existing 
observational evidence reveals that small amplitude oscillations, after excited, are damped 
in short spatial and temporal scales by some as yet not well determined physical mecha- 
nism^). Commonly, these oscillations have been interpreted in terms of linear magnetohy- 
drodynamic (MHD) waves, and this paper reviews the theoretical damping mechanisms that 
have been recently put forward in order to explain the observed attenuation scales. These 
mechanisms include thermal effects, through non-adiabatic processes, mass flows, resonant 
damping in non-uniform media, and partial ionization effects. The relevance of each mech- 
anism is assessed by comparing the spatial and time scales produced by each of them with 
those obtained from observations. Also, the application of the latest theoretical results to 
perform prominence seismology is discussed, aiming to determine physical parameters in 
prominence plasmas that are difficult to measure by direct means. 

Keywords Waves ■ MHD ■ Sun: prominences ■ Sun: filaments • Sun: atmosphere • Sun: 
corona ■ Sun: oscillations 



1 Introduction 

Quiescent solar filaments are clouds of cool and dense plasma suspended against gravity by 
forces though to be of magnetic origin. They form along the inversion polarity line in or 
between the weak remnants of active regions. Early filament observations already revealed 
that their fine s tructure is apparently composed by many horizontal and thin dark threads 
(Engvd3 [l998h . More recent high- resolution H g o bservations obtained with the Swedish 
Solar Telescope (SST) in La Palma fcin et alJ [2005) and the Dutch Open Telescope (DOT) 
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in Tenerife dHeinzel and Anzerll2006t) have allowed to observe this fine structure with much 
greater detail (see llinl 20 id . for a review). The measured average width of resolved thin 
threads is about 0.3 arc. sec (~ 210 km) while their length is between 5 and 40 arc. sec 
(~ 3500 - 28 000 km). The fi ne threads of solar filaments seem to be partially filled with 
cold plasma dlin et al.ll2005h . typically two orders of magnitude denser and cooler than 
the su rrounding corona, and it is generally assumed that they outline their magnetic flux 
tubes tengvoldfl99llLiij2004[Lin et alj200llEngvoldl2008l : lMartin et al.l2008hllln et alj 
120081) . This idea is strongly supported by observations which suggest that they are inclined 
with respe ct to the filament long axis in a similar way to what has been f ound for the mag- 
netic field dLerovl[l9 80; Bom mier et alJfl994l : iBomrnier and Lerovl[l998T) . On the opposite, 
iHeinzel and Anzerl d2006h suggest that these dark horizontal filament fibrils are a projection 
effect. According to this view, many magnetic field dips of rather small vertical extension, 
but filled with cool plasma, would be aligned in the vertical direction and the projection 
against the disk would produce the impression of a horizontal thread. 

Oscillations in prominences and filaments are a commonly observed phenomenon. They 
are us ually classified, in terms o f their amplitude, as small and large amplitude oscilla- 
tions dOliver and Ballesteri 2002 ) . This paper is concerned with s mall amplitude oscilla- 
tions. Large amplitude oscillations have been recently reviewed bv iTripathi eTal] d2009h . 
It is well established that these small amplitude periodic changes are of local nature. The 
detected peak velocity ranges from the noise level (down to 0. 1 km s~ 1 in some cases) to 2- 
3 km s -1 , although larger va lues have also been reported dBashkirtsev and Mashnichll 19841: 



Molowny-Horas et al. 1999). Two-dimensional observations of filaments dYi and Eng vold 



199lUYiet ali i 991) revea ed that individual fibrils or groups of fibrils may oscillate inde 



pendently with their own periods, which range between 3 and 20 minutes. More recently, 
Lin ( 2004) reports spatially coherent oscillations found over slices of a polar crown filament 
covering an area of 1.4 x 54 arc. sec 2 with, among ot her, a significant periodicity at 26 min- 
utes, strongly damped after 4 periods. Furthermore. iLin et al.l d2007h have shown evidence 
about traveling waves along a number of filament threads with an average phase velocity 
of 12 km s _1 , a wavelength of 4" (~ 2800 km), and oscillatory periods of the individ- 
ual threads that vary from 3 to 9 minutes. Oscill a tory event s have been reported from both 
ground-based observations (iTerradas et al]l2002l : Ei^ 2004) , as well as observations from 
instruments onboa rd space-craf ts, such as SoHO felanco et alJI 19991: iRegnier et al1l200ll : 
Pou get et al.ll2006l) and Hinode dOkamoto et al.|[2007l : lNing et al.ll2009l) . The observed peri- 
odic signals are mainly detected from Doppler velocity measure ments and can t herefore be 
associated to the transverse displacement of the fine structures iLin et alJl2009h. Extensive 
review s on small amplit u de oscil l ations i n prominences ca n be found in| 01iver and Ballester 
(20021 ) : lEngvoldl d2004l): Iwiehri d2004l ): IBallesteri (f2006h: iBaneriee et al.l d2007l) : lEngvoldl 
2008 ) : lOliverl d2009r) : IBallesteri feO 1 Ch . and lMackav et alJteOlCh 



Small amplitude oscillations in qui escent filaments have been interpreted in t erms of 
magnetohydrodynamic (MHD) waves loi iver and Ballester 2002; Ballester 2006), which 
has allowed to develop the theoretical models (see IBallesteri 120051 . 120061 . for recent re- 
views). Models of filament threads usually represent them as pa rt of a larger magnetic flux 
tube, whose foot-p oints are anchored in the solar photosphere dBallester and PriestH l989: 
iRempel et al.| [T999). Early works studying filament thread oscillations have considered the 
MHD eigenmodes supported by a filament thread modeled as a Cartesian slab, partially 
filled with p rominence plasma, and embedded in the corona djoarder et al.ll 19971 : iDfaz et al.l 
1200 lL 12003m . Thes e works wer e later extended by considering a more representative cylin- 
drical geometry bv lDfaz et al.l d2002l) . These authors found that the fundamental transverse 
fast mode is always confined in the dense part of the flux tube, hence, for an oscillating 
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cylindrical filament thread, it should be difficult to induce oscillations in adjacent threads, 
unless they are very clos e. Groups of mu lt ithread structures and the ir collective oscillations 
have b een modeled by Diaz et al.l ( 120051) ; iDfaz and Roberts! ( 120061) in Cartesian geometry 
and bv lSoler et al.l ( 12009c ) in cylindrical geometry. 

Time and spatial damping is a recurrently observed characteristic of small amplitude 
prominence oscillations. Observational evidence for the da mping of small amplitude os 



dilations in prominenc e s can be found in [Landman et al. (1977); Tsubaki and Takeuchi 
19861) : iTsubakil Jl988h : Iwiehr et al.l Jl989h: iMolownv-Horas et al] (tl999h; ITerradas etal 



2002 ). and more recently in lLinl J2004h : iBerger et alJ ( l2008h : lNing et all ( 120091) : iLin et al 



20091) . These observational studies have allowed to obtain some characteristic spatial and 



time scales. Reliable val ues for the damping time hav e bee n derived, from differe nt Doppler 
velocity time s erie s, bvlMolownv-Horas et al] d 19991) and lTerradas et all J2002h . in promi- 
nences, and by Lin (2004) in filaments. The values thus obtained are usually between 1 and 4 
times the corresponding period, and large regions of prominences/filaments display similar 
damping times. Se veral th e oretic al mechanisms have been proposed in order to explain the 
observed damping. I Ball ail (2003) estimated, through order of magnitude calculations, that 
several isotropic and anisotropic dissipative mechanisms, such as viscosity, magnetic diffu- 
sivity, radiative losses and thermal conduction cannot in general ex plain the observed wav e 
damping. Linear non - adiab at ic MHD waves have bee n studied by ICarbonell et al] ([2004 ); 
ITerradas et ail feOOll . 120051) : ISoler et al.l J2007L l2008h . The overall conclusion from these 
studies is that thermal mechanisms can only account for the damping of slow waves in 
an efficient manner, while fast waves remain almost undamped. Since prominences can 
be considered as partially ionized plasmas, a possible mechani sm to damp fast waves (a s 
well as Alfven waves) could come from ion-neutral collisions dForteza et al1l2007L [2008), 
although the ratio of the damping time to the period does not completely match the observa- 
tions. Besides non-ideal mechanisms, another possibility t o attenuate fast waves in thin fila- 
ment threads comes from resonant wave damping (see e.g. lGoossens et al" ]|2010l). This phe- 
nomenon is we ll studied for transverse kink waves in coronal loops jGoossens et"aTi r2006 ; 
lGoossens»2 008) and provides a plaus ible explanation for quickly damped transverse lo op os- 
cillations first observed by TRACE ( Aschwanden et al. 19991; iNakariakov et al.lll999l) . The 
time scales of damping produced by these different mechanisms should be compared with 
those obtained from observations. The theoretical approach of many works studying the 
damping of prominence oscillations has been to first study a given damping mechanism 
in simplified uniform and unbounded media thereafter introducing structuring and non- 
uniformity. This has led to an increasing complexity of theoretical models in such a way 
that some of them now combine different damping mechanisms. 

This paper presents results from recent theoretical studies on small amplitude oscilla- 
tions in prominences and the i r dam ping. This top i c has been the theme of several recent 
review works, see e.g. lOlived d2009h and lBallested §010), a fact that reflects the liveliness 
of the subject. Here we extend these works by reporting the last studies and mainly focusing 
on the theory of damping mechanisms for oscillations in prominences and its application 
to prominence seismology, a discipline in a rapidly developing stage, as better observations 
and more accurate theoretical models become available. 



2 Damping of oscillations by thermal mechanisms 



In a seminal paper. lFieldldl965t) studied the thermal instability of a dilute gas in mechanical 
and thermal equilibrium. Using this approach, the time and spatial damping of magneto- 
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hydro dynamic waves in unboun ded media llCarbonell et all [2004. 2009) and in bounded 
slabs dTerradas et alJl200ll 120051) with physical properties akin to those of quies cent solar 
prom i nences , as well as in slabs having prominence-corona physical properties dSoler et all 
l200l [2009a) have been studied. Furthermore, the behavior of non-adiabatic waves in an 
isolated and flowing prominence fine structure dSoler et all2008h has also been considered. 
All these investigati ons have been already reviewed in lOliver! te009h . ballester! told) and 
lMackavetai] ( l2010h . For this reason, in the following we concentrate in a recent investiga- 
tion about the time damping of linear non-adiabatic MHD waves in a multi-thread promi- 
nence when mass flows are also present. 

2.1 Time damping of non-adiabatic magnetoacoustic waves in filament threads with mass 
flows 

In an attempt to explai n the observations about in-phase oscillations of large areas of a fil- 
ament reported bv lLinl ( 12004 ). Ipfaz et alJ d2005h studied the collective oscillations of a set 
of inhomogeneous threads in Cartesian geometry. The configuration was made of five un- 
equally separated threads having different widths and densities, while the magnetic field 
keeps the same strength everywhere. When realistic separations between threads are con- 
sidered, the system oscillates with the only non-leaky mode, and oscillations are in phase 
with similar amplitudes and with a frequency smaller than that of the densest thread. Al- 
though these results show some agreement with observations, the use of Cartesian geometry 
favors collective oscillations since the transversal velocity perturbation has a v ery long tai l 
which allows for an easy interaction between neighboring threads. However, iLinl j2004h 
also pointed out that, after four periods, the oscillations were strongly damped in time and 
that simultaneously flowi ng and osc i llatory structures were also present in filament threads. 
To model this situation, ISoler et all d2009d) have used the T-matrix theory of scattering to 
study the propagation of non-adiabatic magnetoacoustic waves in an arbitrary system of 
cylindrical threads when material flows inside them are present. In this study, the effects of 
radiative losses, thermal conduction and heating have been taken into account, mass flows 
along magnetic field lines have been allowed, and the general case j6 7^ has been consid- 
ered. The equilibrium system is made of two homogeneous and unlimited parallel cylinders, 
prominence threads, embedded in an also homogeneous and unbounded coronal medium 
(see Fig.[T}. For the study, two different configurations have been chosen. 

2.1.1 Identical threads 

In this case the same typical values of prominence density and temperature have been used 
for both threads, as well as typical coronal values for the external density and temperature. 
Furthermore, the magnetic field strengt h is the same e verywhere. In the absence of flow, 
four fast kink modes, already found bv lLunaet all 12008) when studying the collective os- 
cillations of two coronal loops, are present. 

These modes are named as S x ,A x ,S y ,A y , where S and A denote symmetry or antisymme- 
try of the velocity field inside the tubes with respect to the (y, z) plane and the subscripts refer 
to the main direction of polarization of motions. In addition to the kink modes, two further 
fundamental collective wave modes, one symmetric (S z ) and one antisymmetric (A,), mainly 
polarized along the z-direction and corresponding to slow modes have been identified. 

Figure [2^ displays the ratio of the real p art of the frequenc y of the four kink solutions to 
the frequency of the individual kink mode JSoler et alj|2008h . (Ok , as a function of the dis- 
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Fig. 1 Scheme in the xy-plane of the model considered in Sect. 12. l.TI The z- axis is perpendicular to the plane 
of the figure and points towards the reader. Adapted from Soler et al. (20 09d) . 
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Fig. 2 a) Ratio of the real part of frequency, £0r, of the S x (solid line), A x (dotted line), S y (triangles), and 
Ay (diamonds) wave modes to the frequency of the individual kink mode, Oj, as a function of the distance 
between centers, b) Ratio of the damping time to the period versus the distance between centers. Linestyles 
are the same as in panel a). Adapted from Soler et al. 1 2 009d) . 



tance between the center of cylinders, d. It can be seen that the smaller the distance between 
centers, the larger the interaction between threads and the separation between frequencies. 
The frequency of collective kink modes is almost identical to the individual kink frequency 
for a distance between threads larger than 6 or 7 radii. Therefore, we expect the collective 
behavior of oscillations to be stronger when the threads are closer. On the opposite, for larger 
distances, the interaction between threads is much weaker and individual oscillations are ex- 
pected. Furthermore, Figure [2J5 shows the ratio of the damping time to the period, Td/P, of 
the four kink modes as a function of d. This ratio appears to be very large suggesting that 
dissipation by non-adiabatic mechanisms is not very efficient to damp the fast kink modes 
which could be responsible of observed filament threads oscillations. 
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Fig. 3 a) Ratio of the real part of frequency, ft>R, of the S z (solid line) and A z (dotted line) wave modes to 
the frequency of the individual slow mode, co s , as a function of the distance between centers, b) Ratio of 
the damping t ime to the period v ersus the distance between centers. Linestyles are the same as in panel a). 
Adapted from lSoler et~ai] J2009d) . 
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In the case of slow modes, Figure [3^ displays the ratio of the real part of the frequency 
of the S z and A z solutions to the frequency of the individual slow mode, co s . It can be seen 
that the frequencies of the 5- and A z modes are almost identical to the individual slow mode 
frequency, and the strength of the interaction is almost independent of the distance between 
cylinders. This is in agreement with the fact that transverse motions (responsible for the 
interaction between threads) are not significant for slow-like modes in comparison with 
their longitudinal motions. Therefore, the S z and A z modes essentially behave as individual 
slow modes, contrary to kink modes, which display a more significant collective behavior. 
Finally, Figure [3J3 shows Tp/P corresponding to the S z and A z solutions versus d. Slow 
modes are efficiently attenuated by non-adiabatic mecha nisms, with Td/P « 5, which is in 
agreement with previous studies dSoler et al.ll2007L f2008) and consistent with observations. 

Next, the effect of flows on the behavior of collective modes has been studied. Arbitrary 
flows have been assumed in both cylinders while coronal flows have been neglected. First of 
all, we concentrate on transverse modes. We denote the flow in the first cylinder by U\, set- 
ting its value to 20 km s , and we study the behavior of the oscillatory frequency when the 
flow in the second cylinder, denoted by Uj, varies (see Fig.ET l . Sinc e frequencies are almost 
degenerate, we follow the notation of IVan Doorsselaere et al.l J"2008h and call low-frequency 
modes the S x and A v solutions, while high-frequency modes refer toA x and S y solutions. We 
have restricted ourselves to parallel propagation although the argumentation can be easily 
extended to anti-parallel waves. In order to understand the asymptotic behavior of frequen- 
cies in Figure[4] we define the following Doppler-shifted individual kink frequencies: 

Q kl =co k +Uik z , (1) 

Q k2 = co k + U 2 k z , (2) 

Where the + signs arise because of the performed Fourier analysis of the form exp(— ik z z + 
icot). Since U\ is fixed, Q k \ is a horizontal line in Figure[4] whereas Q k 2 is linear with U2- 
Three interesting situations are worth to study and are denoted as a, b, and c in Figure [4] 
(a) When U2 = — 10 km s~', which corresponds to a situation with counter-streaming 
flows, the solutions do not interact with each other and low-frequency (high-frequency) 
solutions are related to individual oscillations of the second (first) thread. For an external 
observer this situation would correspond to an individual thread oscillation. 
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Fig. 4 Ratio of the real part of the frequency, £0r, to the individual kink frequency, fflj., as a function of I] 2 for 
Ui = 20 km s _1 . The solid line corresponds to parallel low-frequency modes (JS X and A v ) while the dashed 
line corresponds to parallel high-frequency solutions (A x and S v ). Dotted lines correspond to the Doppler- 
shifted individual kink frequencies of the threads, Q^i and £2 k 2- The sma ll letters next to the solid line refer 
to particular situations studied in the text. Adapted from Soler et al. 1 2009c). 

(b) When U2 = 20 km s , both flow velocities and their directions are equal in both 
threads. In such a situation, there is a coupling between low- and high-frequency modes, 
and an avoided crossing of the solid and dashed lines is seen in Figure [4] and collective 
oscillations appear. 

(c) When U2 = 27 km s~ ] . This case is just the opposite to (a). Therefore, it corresponds 
again to an individual thread oscillation. 

Considering now slow modes, the behavior of the S z and A z solutions can only be con- 
sidered collective when the flow velocity is the same in both threads because, in such a case, 
both modes couple. When different flows are considered, the S z and A- slow modes behave 
like individual slow modes. The coupling between slow modes is very sensitive to the flow 
velocities, and the S z and A, solutions quickly decouple if U\ and U2 slightly differ. 

2.7.2 Non-identical threads 

Consider, now, a system of two nonidentical threads and focus first on kink modes. From 
above, we expect collective kink motions to occu r when the Doppler- shifted individual kink 
frequencies of both threads coincide. Following I Soler et al . (2009c), the relation between 
flow velocities U\ and U2 for which the coupling takes place is, 

U l -U 2 ^±V2{v A2 -v M ), (3) 

where the + sign is for parallel waves and the — sign is for anti-parallel propagation. A 
similar analysis can be performed for slow modes to obtain, 



U l -U 2 ^±{c s2 -c sl ). 



(4) 



8 



which points out that, in general, the coupling between slow modes occurs at different flow 
velocities than the coupling between kink modes. Therefore, the simultaneous existence of 
collective slow and kink solutions in systems of non-identical threads is difficult. 

Then, we could conclude that the relation between the individual Alfven (sound) speed 
of threads is determinant for the collective or individual behavior of kink (slow) modes. 
In the absence of flows and when the Alfven speeds of threads are similar, kink modes 
are of collective type. On the contrary, when Alfven speeds differ each thread oscillates 
independently. The same happens for slow modes, but replacing the Alfven speeds by the 
sound speeds of threads. 

In summary, when flows are present in the equilibrium, collective motions can be found 
even in systems of non-identical threads by considering appropriate flow velocities. These 
velocities are within the observed values if threads with not too different temperatures and 
densities are assumed. However, since the flow velocities required for collective oscilla- 
tions must take very particular values, such a special situation may rarely occur in real 
prominences. This conclusion has important repercussions for future prominence seismo- 
logical application s, in th e sense that if collective oscillations are observed in large areas of 
a prominence tLinll2004l) . threads in such regions should possess very similar temperatures, 
densities, and magnetic field strengths. 

This study has also confirmed that collective slow modes are efficiently damped by ther- 
mal mechanisms, with damping ratios similar to those reported in observations, Td/P ~ 5, 
while collective fast waves are poorly damped. This is a key point since if the observed ef- 
ficiently damped oscillations are transverse, this could suggest that other mechanisms could 
be at work (see Sects. [3land l4.lt while if the observed oscillations correspond to slow modes, 
then, thermal mechanisms could account for the observed damping. Therefore, it becomes 
crucial to be able to discriminate between the transversal or longitudinal character of the os- 
cillations. However, we should also be aware that because of the inhomogeneous nature of 
filament threads coupling between modes could make a proper mode identification difficult. 



3 Damping of oscillations by ion-neutral collisions 

Since the temperature of prominences is typically of the order of 10 4 K, the prominence 
plasma is only partially ionized. The exact ionization degree of promine nces is unknown and 
the re ported ratio of electron density to neutral hydrogen density (see e.g. lPatsourakos and Vi; 
120021) covers about two orders of magnitude (0.1-10). Partial ionization brings the presence 
of neutrals and electrons in addition to ions, thus collisions between the different species are 
possible. Because of the occurrence of collisions between electrons with neutral atoms and 
ions, and more importantly between ions and neutrals, Joule dissipation is enhanced, when 
compared to the fully ionized case. 

The main effects of partial ionization on the properties of MHD waves are manifested 
through a generalized Ohm's law, which adds some extra terms in the resistive magnetic 
inducti on equation, in com parison to the fully ionized case. The induction equation can be 
cast as llSoler et al.l l2009d) 



— -i = Vx (vi xB )-Vx(?jVxB 1 )-|-Vx{tja[(VxBi) xB ] xB } 




(5) 
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with Bi, and V] the perturbed magnetic field and velocity, respectively. Quantities 77, 77 a, 
and rju in Equation (|5} are the coefficients of ohmic, ambipolar, and Hall's magnetic dif- 
fusion, and govern the collisions between the different plasma species. Ohmic diffusion is 
mainly due to electron-ion collisions, ambipolar diffusion is mostly caused by ion-neutral 
collisions, and Hall's effect is enhanced by ion-neutral collisions since they tend to de- 
coupl e ions from the magnetic field while electrons remain able to drift with the magnetic 
field dPandev and Wardl e 2008). The ambipolar diffusivity can be expressed in terms of the 
Cowling's coefficient, Tjc, as 

The quantities 77 and 77c correspond to the magnetic diffusivities longitudinal and perpen- 
dicular to magnetic field lines, respectively. For a fully ionized plasma, 77c = 77 , there is no 
ambipolar diffusion, and magnetic diffusion is isotropic. Due to the presence of neutrals, 
"Hc^ "H, meaning that perpendicular magnetic diffusion is much more efficient than longi- 
tudinal magnetic diffusion in a partially ionized plasma. It is important to note that 77c 2> 77 
even for a small relative density of neutrals. 



3.1 Partial ionization effects in a homogeneous an unbounded prominence medium 
3.1.1 Time damping of magnetohydrodynamic waves 

Several studies have c onsidered the damping of MHD waves in the solar atmosphere in par- 
tially ionized plas mas ( IDe Pontieu et alj|200ll ; |james et al.ll2003Kl Khodache nko et alj 2004; 
iLeake et al.ll2005h . In the context of solar prominences iForteza et al.1 112007") derived the full 
set of MHD equations for a partially ionized, one-fluid hydrogen plasma and applied them 
to the study of the time damping of linear, adiabatic fast and slow magneto-acoustic waves 
in an unbounded prominence medium. A partially ionized plasma can be represented as a 
single-fluid in the strong coupling approximation, which is valid when the ion density in 
the plasma is small and the collision time between neutrals and ions is short compared with 
other timescales in the problem. Using this approximatio n we can describ e the very low 
frequency and large -scale fluid-like behaviour of plasmas dGoossensl l2003). The study by 
IForteza et al.l |2007) was later extended to the non- adiabatic case, inclu ding thermal conduc- 
tion by neutrals an d electrons and radiat ive losses ( IForteza et al.f 2008). The most importan t 
res ults obtained by 



Forte za et al.l 12 007) have been summarized recently by lOlive j |2009) 



and iBallesteil (12010). IForteza et al.l ( 12007 ) consider a uniform and unbounded prominence 



plasma and find that ion-neutral collisions are more important for fast waves, for which the 
ratio of the damping time to the period is in the range 1 to 10 5 , than for slow waves, for 
which values in between 10 4 and 10 8 are obtained. Fast waves are efficiently damped for 
moderate values of the ionization fraction, while in a nearly fully ionized plasma, the small 
amount of neutrals is insufficient to damp the perturbations. 

In the above studies, a hydrogen plasma has been considered, however, 90% of the 
prominence chemical composition is hydrogen while the remaining 10% is helium. There- 
fore, it is of great interest to know the effect of the presence of helium on the behavior 
of magnetohydrodynamic waves in a p artially ionized pla sma with prominence physical 
propert ies. This study has b een done by ISoler et a l. ( 2009e) in a medium like that consid- 
ered in IForteza et al but composed of hydrogen and helium. The species present 
in the medium are electrons, protons, neutral hydrogen, neutral helium (Hel), and singly 
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ioniz ed helium (Hell), while the presence of He III is negligible dGouttebroze and Labrossd 
2009). Under such conditions the basic MHD equations for a non-adiabatic, partially ion- 
ized, single-fluid plasma have been generalized. 

The hydrogen ionization degree is characterized by p.n which varies between 0.5, for 
fully ionized hydrogen, and 1 for fully neutral hydrogen. The helium ionization degree is 
characterized by 5n e = fj^, where (§Heii and <§Hei denote the relative densities of single 
ionized and neutral helium, respectively. Figure [5] displays To /P as a function of k for the 
Alfven, fast, and slow waves, and the results corresponding to several helium abundances 
are compared for hydrogen and helium ionization degrees of p.n =0.8 and 5n e = 0.1, re- 
spectively. We can observe that the presence of helium has a minor effect on the results. 
In the case of Alfven and fast waves (Figs.[5t,b), a critical wavenumber, k*, at which the 
real part of the frequency becomes zero occurs. In a partially ionized plasma, this critical 
wavenumber, k*, is given by 

p = 2v 'a m 

c cos d(j] c + ri tan 2 0) ' 

where va is the Alfven speed and 9 the angle between the wavevector and the equilibrium 
magnetic field. The Cowling's diffusivity, rjc, depends on the fraction of neutrals and the 
collisional frequencies between electrons, ions and neutrals. When a fully ionized plasma 
is considered the numerical value of both diffusivities is the same and in a fully ionized 
ideal plasma this numerical value is taken equal to zero, therefore, the critical wavenumber 
goes to infinity. When partially ionized plasmas are considered, wavenumbers greater than 
the critical value only produce purely damped perturbations. Since Cowling's diffusivity 
is larger in the presence of helium because of additional collisions of neutral and singly 
ionized helium species, is shifted toward slightly lower values than when only hydrogen 
is considered, so the larger ^Hei, the smaller k%. In the case of the slow wave (Fig.[3J;), the 
maximum and the right hand side minimum of Tn/P are also slightly shifted towa rd lower 
values of k. Previous results from ICarbonell e t al. (200'}) and lForteza et alJ J2008h suggest 
that thermal conduction is responsible for these maxima and minima of To /P- The additional 
contribution of neutral helium atoms to thermal conduction produces this displacement of 
the curve of Tp/P. In the case of A l fven a nd fast waves, this effect is not very important. 
Although Gouttebroze and Labrossd ( 120091) suggest that for central prominence temperature 
a realistic ratio between the number densities of <§Heii to <§e e i is 10%, in Figure [5] and for 
sake of comparison, the results for <§e e i — 10% and 5n e = 0.5 have been also plotted. 

Finally, the thermal mode has been considered. Since it is a purely damped, non-propa- 
gating disturbance (cor = 0), we only plot the damping time, To, as a function of k for flu = 
0.8 and ^He = 0. 1 (Fig.[5Jl). We observe that the effect of helium is different in two ranges of 
k. For k > 10~ 4 m _1 , thermal conduction is the dominant damping mechanism, so the larger 
the amount of helium, the smaller To because of the enhanced thermal conduction by neutral 
helium atoms. On the other hand, radiative losses are more relevant for k < 10~ 4 m _1 . In 
this region, the thermal mode damping time grows as the helium abundance increases. Since 
these variations in the damping time are very small, we again conclude that the damping 
time obtained in the absence of helium does not significantly change when helium is taken 
into account. 

In summary, the study bv lSoler"et ail d2009eh points out that the consideration of neu- 
tral or single ionized helium in partially ionized prominence plasmas does n ot modify the 
behav ior o f linear, adiabatic or non-adiabatic MHD waves already found by iForteza et alJ 
d2007l) and lForteza etafl d2008l) . 



11 




fe (rrC ) k (m ' ) 

Fig. 5 (a)-(c): Ratio of the damping time to the period, To/P, versus the wavenumber, k, corresponding to 
the Alfven wave, fast wave, and slow wave, respectively, (d) Damping time, Tj, of the thermal wave versus 
the wavenumber, k. The different linestyles represent: ^Hel — 0% (solid line), i^Hcl = 10% (dotted line), and 
^Hel = 20% (dashed line). In all computations, /2h = 0.8, 8u e = 0.1, and the angle between the wavevector 
and the equilibrium magnetic field is 6 = Jt/4. The results for i^Hcl = 10% and <5h c = 0.5 are plotted by means 
of symbols for comparison. The shade d regions correspond to the range of typically observed wavelengths of 
prominence oscillations. Adapted from Soler et al. 1 20 09a) . 



3.1.2 Spatial damping of magnetohydrodynamic waves 

iTerradas et aT] £2002) analyzed small amplitude oscillations in a polar crown prominence 
and reported the presence of a plane propagating wave as well as a standing wave. In the 
case of the propagating wave, which was interpreted as a slow MHD wave, the amplitude of 
the oscillations spatially decreased in a substantial way after a distance of 2 — 5 x 10 4 km 
from the location where wave motion was being generated. This distance can be considered 
as a typical spatial damping length, L^, of the oscillations. On the other hand, a typical 
feature in promi nence oscillations is the presence of flows which are observed in H a , UV 
and EUV lines ( lLabrosse etal]|2010l). In Hq quiescent filaments , the observed velocities 



range from 5 to 20 km s 1 ( Zirk er et al Jl998l :l iLin et al.l2 003. 2007) and, because of physical 



conditions in filamen t plasma, they seem to be field-aligned. Recently, observations made 
with Hinode/SOT by Oka moto et al] d2007h reported the presence of synchronous vertical 
oscillatory motions in the threads of an active region prominence, together with the presence 
of flows along the same threads. However, in limb prominences differen t kinds of flows are 
observed and, for instance, observations made by Berger et al ] d2008h with Hinode/SOT 
have revealed a complex dynamics with vertical downflows and upflows. 

The spatial damping of magneto hydrodynamic w a ves in a homogeneous and unbounded 
fully ionized plasma was studied bv lCarbonell et alJ 120061) . Recently, the spatial damping 
of linear non-adiabatic magnetohydrodynamic waves in a ho mogenous, unbounded, magne- 
tized and flowing partially ionized plasma has been studied bv lCarbonell et al.l d2010h . Since 
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Fig. 6 Damping length, wavelength, and ratio of 
the damping length to the wavelength versus pe- 
riod for the three (solid, dashed, dotted) Alfven 
waves in a partially ionized plasma with an ion- 
ization degree fl — 0.8 and with a background 
flow of 10 km s . In all the panels, the shaded 
region corresponds to the interval of observed pe- 
riods in prominence oscillations. 



we consider a medium with physical properties akin to those of a solar prominence, the den- 
sity is po = 5 x KT 11 kg m~ 3 , the temperature 7b = 8000 K, the magnetic field |Bq|= 10 
G, and, in general, a field-aligned flow with vo = 10 km s~' has been considered. 
The dispersion relation for Alfven waves with a background flow is given by: 

»'o(f]ccos 2 6 + r\ sin 2 0)& 3 
+ [( v o ~ v a~ !r ?c©) cos 6 — ir\w sin 6 tan 6] k 2 — Icovok 

+ffl 2 sece = 0, (8) 

with d the angle between the wavevector and the equilibrium magnetic field. The increase 
in the degree of the dispersion relation, with respect to the typical dispersion relation for 
Alfven waves, is produced by the joint presence of flow and resistivities. In this case, we 
should obtain three propagating Alfven waves, and Figure [6] shows the numerical solution 
of dispersion relation ((8). For the entire interval of periods considered, a strongly damped 
additional Alfven wave appears, while on the contrary, the other two Alfven waves are very 
efficiently damped for periods below 1 s. However, within the interval of periods typically 
observed in prominence oscillations these waves are only efficiently attenuated when almost 
neutral plasmas are considered. 

The dispersion relation for thermal and magnetoacoustic waves in presence of a back- 
ground flow is given by 

(n 2 -k 2 A 2 ){ik 2 rj C n-n 2 )+k 2 vl{n 2 -k 2 f A 2 )+ik 2 k 2 v 2 A A 2 Ep n =o, (9) 

where A 2 is the non-adiabatic sound speed squared dForteza et al.ll2008l : ISoler et al. 2008) 
and is defined as 
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Fig. 7 Damping length, wavelength, and ratio of the damping length to the wavelength versus period for the 
non-adiabatic fast (left panels), slow (right panels) waves in a partially ionized plasma with an ionization 
degree ft = 0.8 (solid) and fl = 0.95 (dashed). The flow speed is 10 km s _1 , 



where A and H, including optically thin radiative losses, thermal co nduction by electron s 
and neutrals, and a constant heating per unit volume, were defined in lForteza et alj ([2008). 
When Equation $9} is expanded, it becomes a seventh degree polynomial in the wavenumber 
k, whose solutions are three propagating fast waves, two slow waves and two thermal waves. 

Figure[7]displays the behavior of the damping length, wavelength and ratio of damping 
length versus wavelength for two of the three fast waves and slow waves. The curves for the 
third fast wave, which owes its existence to the joint action of flow and resistivity are similar 
to those corresponding to the strongly damped Alfven wave in Figure [6] The most interest- 
ing results are those related with the ratio L&/X. For fast waves, this ratio decreases with 
the period becoming small for periods below 10~ 2 s while for one of the slow waves, the 
ratio becomes very small for periods typically observed in prominence oscillations. When 
ionization is decreased, slight changes of the above described behavior occur, the most im- 
portant being the displacement towards longer periods of the peak of most efficient damping 
corresponding to slow waves. 

In summary, Alfven waves in a partially ionized plasma can be spatially damped. When 
the ionization decreases, the damping length of these waves also decreases and the efficiency 
of their spatial damping in the range of periods of interest is improved, although the most 
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efficient damping is attained for periods below 1 s. A new feature is that when a flow is 
present a new third Alfven wave, strongly attenuated, appears. The presence of this wave 
depends on the joint action of flow and resistivities, and it could only be detected by an 
external observer to the flow. 

In the case of non-adiabatic magnetoacoustic waves, when partial ionization is present 
the behavior of fast, slow and thermal waves is strongly modified. In particular, the damping 
length of a fast wave in a partially ionized plasma is strongly diminished by neutrals thermal 
conduction for periods between 0.01 and 100 s, and, at the same time, the radiative plateau 
present in fully ionized ideal plasmas almost disappears. The behavior of slow waves is not 
so strongly modified as for fast waves, although thermal conduction by neutrals also dimin- 
ishes the damping length for periods below 10 s, and a short radiative plateau still remains 
for periods between 10 and 1000 s. Thermal waves are only slightly modified although the 
effect of partial ionization is to increase the damping length of these waves, just the oppo- 
site to what happens with the other waves. Next, when a background flow is included, a new 
third fast wave appears which, again, is due to the joint action of flow and resistivities. Also, 
in the presence of flow, wavelengths and damping lengths are modified, and since for slow 
waves sound speed and observed flow speeds are comparable this means that the change in 
wavelength and damping length are important, leading to an improvement in the efficiency 
of the damping. Moreover, the maximum of efficiency is displaced towards long periods 
when the ionization decreases, and for ionization fractions from 0.8 to 0.95 it is clearly lo- 
cated within the range of periods typically observed in prominence oscillations with a value 
of Ld/A smaller than 1. This means that for a typical period of 10 3 s, the damping length 
would be between 10 2 and 10 3 km, the wavelength around 10 3 km and, as a consequence, 
in a distance smaller than a wavelength the slow wave would be strongly attenuated. 

In conclusion, the joint effect of non-adiabaticity, flows and partial ionization allows to 
damp slow waves in an efficient way within the interval of periods typically observed in 
prominences. 

3.2 Partial ionization effects in a cylindrical filament thread model 

Recently, ISoler et all d2009bh have applied the equations derived by iForteza et al.l l l2007h 
to the study of MHD waves and their time damping in a partially ionized filament thread. 
The adopted thread model is the classic one-dimensional magnetic flux tube, with promi- 
nence conditions, embedded in an unbounded medium with coronal conditio ns. A uniform 
magne tic field, Bq pointing along the axis of the tube is considered. As in IForteza et al.l 



(2007), the one fluid approximation for a hydrogen plasma is considered. The internal and 
external media are characterized by their densities, temperatures, and the relative densities 
of neutrals, ions, and electrons. The contribution of the latter is neglected. The ionization 
fraction is thus defined as p = 1/(1 + Q), with the relative density of ions. For a fully 
ionized plasma p = 0.5 (£,- = 1), while for a neutral plasma p = 1 (£,- = 0). The external 
coronal medium is considered as fully ionized, while the ionization fraction in the internal 
filament plasma is allowed to vary in this range. Note that any value of p outside this range 

is physically meanin gless. 

In their analysis ISoler et al. (2009b) neglect Hall's term, which is important only for 
frequencies larger than ~ 10 4 Hz, much larger than the observed frequencies of promi- 
nence osci llations (for a dimensional analysis that further justifies this approximation see 
ISoler et aTll2009dh . After Fourier analyzing the linear MHD wave equations by assuming 
perturbations of the form exp(iwt + im(p — ;'fc,z). ISoler et al.l d2009bl) note that terms with rjc 
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appear accompanied by longitudinal derivatives, while terms with r\ correspond to radial and 
azimuthal derivatives. By defining and as typical length-scales parallel and perpen- 
dicular to the magnetic field, the first is associated to the wavelength of perturbations in the 
longitudinal direction, X z ~ 2nk z , while the second is related to the radius of the structure, 
a. This allows to define the corresponding Reynolds numbers in the parallel and perpendic- 
ular directions as R m \\ = c s fa/r\ and R m ± = 47t 2 c s f/r)ck z a, where the typical velocity scale 
has been associated to the sound speed in the filament, c s j. The parallel Reynolds number 
is independent of the wavenumber, while the relative importance of Cowling's diffusion in- 
creases with k z . A simple calculation reveals that Cowling diffusion is dominant for the cases 
of interest. Consider, for instance, k z a = 1, a = 100 km, and p,f = 0.8. Then R m ii = 7 • 10 6 
and R m ± =4 - 1 2 . The wavenumber for which both ohmic and Cowling's diffusion have the 
same importance can be estimated to be X ~ [5 • 10 3 — 10 5 ] km, for a 6 [75 — 375] km. In the 
range of observed wavelength s (k 7 a ~ [10~ 3 — 1 1 ]) both Cowling's and ohmic diffusion 
could therefore be important. Sole r et al.l d2009bl) analyze separately the effect(s) of partial 
ionization in Alfven, fast kink, and slow waves. 

Alfven waves are incompressible disturbances with velocity perturbations polarized in 
the azimuthal direction. In a resistive medium these velocity pertur bations are not strictly 
confined to magnetic surf aces, but have a global nature (see e.g., iFerraro and PlumptorJ 
Il96ll) . ISoler et alJ |2009b) show that, by considering solutions with m = (no azimuthal 
dependence) and defining a modified Alfven speed squared as = v\ + iwrjc, with the 
ideal Alfven speed, the azimuthal components of the momentum and induction equations 
can be combined to obtain a Bessel equation of order one for the perturbed magnetic field 
component in the azimuthal direction. By solving this equation, the phase speed and damp- 
ing rate of Alfven waves can be studied as a function of the wavelength of perturbations 
for different ionization degrees and values of the ohmic dissipation. It turns out that Alfven 
wave propagation is constrained between two critical wavenumber values. These critical 
wavenumbers are, however, outside the range that corresponds to the observed wavelengths 
(Fig.[8]top). The small critical wavenumber is found to be insensitive to the ionization frac- 
tion, while the large critical wavenumber strongly depends on this parameter. The obtained 
values of the damping time over the period are independent of the ionization degree, for 
small wavenumber values, while they are affected for large wavenumber values. They are 
found to be in between 10 and 10 2 times the corresponding period, in the range of typically 
observed wavelengths. By considering solutions to the dispersion relation by neglecting 
separately one of the two poss ible damping mechanisms, i.e., partial ionization and ohmic 
dissipation, ISoler et al .(2009b) observe that the presence of neutrals has an important effect 
on the damping time for large wavenumbers, while ohmic diffusion dominates for small 
wavenumbers (Fig. [9^). 

The propagation of transverse kink waves is also found to be constrained by two critical 
wavenumbers (Fig. [8] middle), a result that can be easily understood by noting that both in 
the short- and long-wavelength regimes the asymptotic approximations for the kink mode 
frequency directly depend on the Alfven speed. Within the relevant range of observed wave- 
lengths, the phase speed closely corresponds to the ideal counterpart, = Co/k z , so non- 
ideal effects are irrelevant to wave propagation for the observed wavelengths. The behavior 
of the damping rate as a function of wavelength and ionization fraction is seen to closely re- 
semble the result obtained for Alfven waves, with values of ij jP above observed values, in 
the range of observed wavelengths. Therefore, neither ohmic diffusion nor ion-neutral col- 
lisions seem to provide damping times as short as those observed in transverse kink waves 
in filament threads. Only for an almost neutral plasma, with p-f > 0.95 and very short wave- 
lengths the obtained damping rates are compatible with the observed time-scales. As for 
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Fig. 8 Phase speed (left panels) in units of the Alfven speed, kink speed and internal cusp speed from top 
to bottom, and ratio of the damping time to the period (right panels) as a function of k z a for Alfven (top 
panels), kink (middle panels) and slow (bottom panels) waves. In all panels different linestyles represent 
different ionization degrees: U/ = 0.5 (dotted), flf = 0.6 (dashed), flf = 0.8 (solid), and flf = 0.95 (dash- 
dotted). Symbols are the approximate solution given by Equation (36) in lSoler et all <2009bl) for fl f = 0.8. 
The shaded z ones correspond to the range of typically observed wavelengths of prominence oscillations. 
Adapted from lSoler et~aTM2009bl) . 



Alfven waves, similar results are obtained for the kink mode when comparing the damp- 
ing rate by neglecting ohmic diffusion or ion-neutral collisions separately (Fig.|9p). Ohmic 
diffusion dominates for small wavenumbers, while ion-neutral collisions are the dominant 
damping mechanism for large values of k z a. 

Finally, ISoler et alj j2009bh analyze the propagation properties and the damping by ion- 
neutral collisions for slow waves. The analysis concentrates on the radially fundamental 
mode with m = 1, since the behavior of the slow mode is weakly affected by the value 
of the azimuthal wavenumber. It turns out that slow wave propagation is constrained by 
one critical wavenumber, which strongly depends on the ionization fraction, in such a way 
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Fig. 9 Ratio of the damping time to the period as 
a function of k z a corresponding to (a) the Alfven 
mode, (b), the kink mode, and (c) the slow mode, 
with fif = 0.8. Solid lines are the complete solu- 
tion considering all terms in the induction equa- 
tion. The dotted and dashed lines are the results 
obtained by neglecting ohmic diffusion (T) = 0) 
or ion-neutral collisions (rjc = rj), respectively. 
The vertical dot-dashed lines indicate the analyt- 
ically estimated transitional wavenumber values 
betwe en both regimes of dominance. Adapted 
from lSoler etafl |2009bl) . 



that for k- < k zcr \x the wave is totally damped. More importantly, for large enough values 
of p., the corresponding critical wavelength lies in the range of observed wavelengths of 
filament oscillations (Fig.[8]bottom). As a consequence, the slow wave might not propagate 
in realistic thin filament threads. By computing the damping rate, it is found that ion-neutral 
collisions are a relevant damping mechanism for slow waves, since very short damping 
times are obtained, especially close to the critical wavenumber. By comparing the separate 
contributions of ohmic diffusion and ion-neutral collisions, the slow mode damping is seen 
to be completely dominated by ion-neutral collisions (Fig. |9};). Ohmic diffusion is found to 
be irrelevant, since the presence of the critical wavenumber prevents slow wave propagation 
for small wavenumbers, where ohmic diffusion would start to dominate. 



4 Resonant damping of filament thread oscillations 

As discussed in the previous sections, non-adiabatic MHD waves and partial ionization seem 
to be able to explain the time damping of oscillations in the case of slow waves. The question 
remains about what physical mechanism(s) could be responsible for the rapid time damping 
of transverse kink waves in thin filament threads, which seem to be rather insensi tive to ther- 
mal effe cts. Partial ionization could be relevant, in view of the results obtained bv lSoler et alj 
(2009b), but the ratios of the damping time to the period obtained for a cylindrical filament 
thread are still one or two orders of magnitude larger that those observed. The phenomenon 
of resonant wave damping in non-uniform media has provided a pl ausible explanation in 
the context of quickly damped transv erse coronal loop oscillations dGoossens et alJl2~002l : 
iGoos sens 20081; iGoossens et"aT]|2010l) . This mechanism relies in the wave coupling and en- 
ergy transfer between oscillatory modes due to the spatial variation of physical properties 
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that in turn define a non-uniform Alfven and/or cusp speed. Because of the highly inhomo- 
geneous nature of quiescent filament threads at their transverse spatial scales, it is natural to 
believe that resonant damping must be operative whenev er transverse k ink waves propagate 
al ong those struc t ures. This idea has been put forward bv lArregui et al. I l2008h . The analysis 
bv lArregui et"aT] ( l2008h is restricted to the damping of kink oscillations due to the resonant 
coupling to Alfven waves in a pressureless (zero plasma-/?) plasma. It has bee n extended 
to the case in which both the Alfven and the slow resonances are present by ISoler et all 
The main results from these two studies are summarized. 



4. 1 Resonant damping in the Alfven continuum 

Given the relatively simple structure of filament threads, when compared to the full promi- 
nence/filament system, the magnetic and plasma configuration of an indivi dual and iso- 
lated t hread can be theoretically approximated using a rather simplified model. lArreeui et alj 
(2008) make use of a one-dimensional, straight flux tube model in a gravity-free environ- 
ment. In a system of cylindrical coordinates (r, <p, z), with the z-axis coinciding with the axis 
of the tube, the uniform magnetic field is pointing in the z-direction, B = Se z . As gas pres- 
sure is neglected slow modes are absent and the oscillatory properties of the remaining fast 
and Alfven MHD waves and their mutual interaction is analyzed. In particular the analysis 
concentrates on the fundamental kink mode. In such a straight field configuration the zero 
plasma- jS approximation implies that the field strength is uniform and that the density pro- 
file can be chosen arbitrarily. The non-uniform filament thread is then modeled as a density 
enhancement with a one-dimensional non-uniform transverse distribution of density, p(r), 
across the structure. The internal filament plasma, with uniform density, py, occupies the 
full length of the tube and is connected to the coronal medium, with uniform density, p c , by 
means of a non-uniform transitional layer of thickness /. If a denotes the radius of the tube, 
the ratio l/a provides us with a measure of the transverse inhomogeneity length-scale, that 
can vary in between l/a = (homogeneous thread) a nd I /a = 2 (ful l y non -uniform thread). 
The explicit expression for the density profile used by Arregu Tet alJ 1 20081) is 



with 



r < a — 1/2, 

Po0)= { Ptr(r), if a-l/2<r<a + l/2, (11) 
r> a + l/2, 



p t (r) = ¥Ul+^)-fl-^U[5(r- fl )|k (12) 




2 IV Pf) \ Pf) L ' 



Typical values for the filament and coronal densities are py = 5 x 10 11 kg m 3 and p c = 
2.5 x 10~ 13 kg m~ 3 , the density contrast between the filament and coronal plasma being 
Pflp c = 200. 

Observations of transverse oscillations in filament threads can be interpreted in terms 
of linear kink waves. When considering perturbations of the form f(r)exp(i(cot + m(p — 
k z z)), with m and k z the azimuthal and longitudinal wavenumbers and CO the oscillatory 
frequency, the fundamental kink mode has m = 1, and produces the transverse displacement 
of the tube as it propagates along the density enhancement. This mode is therefore consistent 
with the detected Doppl er velocity measure ments and the associated transverse swaying 
motions of the threads dlin et alii 20071 , 12009). In the absence of a non-uniform transitional 
layer, i.e., l/a = 0, the density is uniform in the internal and external regions and changes 
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discontinuously at the tube radius a. Then, a well known dispersion relation is obtained by 
imposing the continuity of the radial displacement, | r , and the total pressure perturbation, 
Pt, at r = a. This dispersion relation (Edwin & Roberts 1983) is 



D m (a,k t ) = J^-J?i = (13) 

r T,e r T,i 



where the indices "i" and "e" refer to internal and external, respectively, and the prime 
denotes a derivative with respect to the radial direction. In the commonly used thin tube 
or long wavelength approximation (k : a << 1), the kink mode frequency can be calculated 
explicitly as 



®k = h\ n ' n » ( 14 ) 

Pf+Pc 

with VAf.c = B/y/fipf :C the filament and coronal Alfven velocities. The period of kink os- 
cillations with a wavelength A = 2k /k z can be written, in terms of the density contrast, 
as 



2 V M 

Note that the factor containing the density contrast varies between v2 and 1, when c is 
allowed to vary between a value slightly larger that 1 (extremely tenuous thread) and c — > °°. 
This has consequences in the application of the model to perform prominence seismology 
(see Sect. 6). 

In Equation ( 114b the eigen-frequency of the fundamental kink mode is in between the 
internal and external Alfven frequencies, C0a/ < < COa c , with <8a/ c = k z VAf. c - This means 
that when the discontinuous jump in density is replaced by a continuous variation in a non- 
uniform layer, of thickness /, going from P/ to p c , the fundamental kink mode has its eigen- 
frequency in the Alfven continuum, and thus couples to an Alfven continuum mode. This 
results in a transfer of wave energy from the transverse motion of global nature to azimuthal 
motions of localized nature, and the time damping of the kink mode. Asymptotic analyt- 
ical expressions for the damping time, Tj, can be obtained under the assumption that the 
transverse inhomogeneity length-scale is small (/ /a <§; 1). This is the so-called thin bound- 
ary approximation, which assumes that the thickness of the dissipative layer, namely 5a, 
coincides with the width of the inhomogeneous transitional layer. This condition is approxi- 
mately verified for thin layers, i.e., l/a <C 1, w hich makes the appro x imation very accurate i n 
such a case. The method w as outlined by e.g. lSakurai et al 1 <199lh : lGoossens et alj fl995b : 
iTirrv and Goossensl \ 19961) and makes use of jump conditions to obtain analytical expres- 
sions for the dispersion relation. For the purposes of our description it will suffice to write 
down these conditions. In a straight and constant magnetic field, we have for the Alfven 
resonance 

2/2 

[| r ] = -in ™ ' A p T , [pr] = 0, at r = r A , (16) 

<t>A\ d rPo\A 

where o>a and |<9 r PoU are tne Alfven frequency and the modulus of the radial derivative of 
the density profile, evaluated at the resonant point. These jump conditions allow us to obtain 
a dispersion relation of the form 
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Fig. 10 (a)-(c): Damping time over period for fast kink waves in filament threads with a = 100 km. In all 
plots solid lines correspond to analytical solutions given by equation (18) , with F = (a): As afunction of 
density contrast, with l/a = 0.2 and for two wavelengths, (b): As afunction of wavelength, with I /a = 0.2, and 
for two density contrasts, (c): As a function of transverse inhomogeneity length-scale, for two combinations 
of wavelength and density contrast, (d) Percentage difference, A, with respect to analytical formula U8t for 
different combinations of wavelengt h, X = 30a (dashed lines); A = 200a (dash-dotted lines), and density 
contrast. Adapted from Arregui et al. 1 2008). 



D(a>,k z ) = -™ m2/ ? A * P \ ■ 07) 

(Oj \d r po\A 

When the long wavelength and the thin boundary approximations are combined, the ana- 
lytical expressi on for the damping time oyer period for the kink mode (m = 1) can be writ- 
ten as (see e.g. Hollweg and Yangl fl988l ; ISakurai et al.ll 199 ll : iGoossens et al.l l 1991 . 1 19951 : 



Ruderman and Robert sll2002h 



Here F is a numerical factor that depends on the particular variation of the density in the 
non-u niform layer. For a linear variation F = 4/ti 2 ([Hollweg and Yandl98ctlGoossens et al.l 
1 1 992h ; for a sinusoidal variation F = 2/ 71 dRuderm an and Robertsl l2002h . We can already 
anticipate that, for example, considering c = 200 as a typical density contrast and I /a = 0.1 
Equation d 1 8b predicts a damping time of ~ 6 times the oscillatory period. 

Figure [10] shows analytical estimates computed by lArreeui et al.l ( 120081) using Equa- 
tion d 1 8b (solid lines). The damping is affected by the density contrast in the low contrast 
regime and Tj/P rapidly decreases for increasing thread density (Fig. 1 10b). Interestingly, 
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it stops being dependent on this parameter in the large contrast regime, typical of filament 
threads. The damping time over period is independent of the wavelength of perturbations 
(Fig. [TOb). but rapidly decreases w ith increasing inhomogeneity length-scale (Fig. [TOb). 
These results led lArregui et al.l d2008h to propose that resonant absorption in the Alfven con- 
tinuum is a very efficient mechanism for the attenuation of fast waves in filament threads, es- 
pecially because large thread densities and transverse plasma inhomogeneities can be com- 
bined together. 

The analysis bv lArreguiet~ai](l2008l) is completed by computing numerical approxima- 
tions to the solutions outside the thin tube and thin boundary approximations used to derive 
Equations dl5b and d 1 8b . which may impose limitations t o the applicability of the obtained 
results to filament thread oscillations. lArregui et al.l ( 120081) find that analytical and numerical 
solutions display the same qualitative behavior with density contrast and transverse inhomo- 
geneity length-scale (Figs.fTOh. c). Now the damping time over period slightly depends on 
the wavelength of perturbations (Fig. 1 10b). Equation dl 8b underestimates/overestimates this 
magnitude for short/long wavelengths. The differences are small, of the order of 3% for 
c = 10, and do not vary much with density contrast for long wavelengths (A = 200a), but in- 
crease until 6% for short ones (A = 30a). The long wavelength approximation is responsible 
for the discrepancies obtained for thin non-uniform layers (Fig. 1 10b). Figure [TOtl shows how 
accurate Equation d 1 3b is for different combinations of wavelength, density contrast, and in- 
homogeneity length-scale. For thin layers (I /a = 0.1) the inaccuracy of the long wavelength 
approximation produces differences up to ~ 10% for the combination of short wavelength 
with high contrast thread. For thick layers, differences of the order of 20% are obtained. 
Here, the combination of large wavelength with high contrast thread produces the largest 
discrepancy. 

Numerical results allow the computation of more accurate values, but do no change the 
overall conclusion regarding the efficiency and properties of resonant damping of transverse 
oscillations in filament threads. Resonant damping in the Alfven continuum appears as a 
very efficient mechanism for the explanation of the observed damping time scales. 



4.2 Resonant damping in the slow continuum 

Although the plas ma jS in solar prominences is probably small, it is definitely nonzero. 
ISoler et al.U2009d) have recen tly incorporated gas pressure to the cylindrical filament thread 



model of lArregui et al. (2008). This introduces the slow mode in addition to the kink mode 
and the Alfven continuum. In the context of coronal loops, which are presumably hotter and 
denser than the surrounding corona, the ordering of sound, Alfven and kink speeds does 
not allow for the simultaneous matching of the kink frequency with both Alfven and slow 
continua. Because of their relatively higher density and lower temperature conditions, this 
becomes possible in the case of filament threads. Therefore, the kink mode phase speed is 
also within the slow (or cusp) continuum, which extends between the internal and exter- 
nal sound speeds, in addition to the Alfven continuum. This brings an addition al damping 
mechan ism. Its contribution to the total resonant damping has been assessed bv lSoler et al.l 
k009d) . 

The study by ISoler et al. (20091) considers a filament thread model equivalent to that 



of lArregui et al.l d2008[) . i.e., a straight cylinder with prominence-like conditions embedded 
in an unbounded corona, with a transverse inhomogeneous layer between both media. The 
same one-dimensional density profile given by Equation i ll U is used. The internal and exter- 
nal densities are set to p,- = 5 x 10~ n kg m~ 3 and p e = 2.5 x 10~ 13 kg m~ 3 , giving a density 
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contrast of pi/p e = 200. The plasma temperature is related to the density through the usual 
ideal gas equation and values of 7} = 8000 K and T e = 10 6 K are taken. The magnetic field, 
pointing in the axial direction is taken to be uniform, with a value So = 5 G everywhere. 
Under this conditions, the plasma-/? ~ 0.04. 

In order to compute the analytical contribution to the damping of the kink mode due to 
the slow resonance a similar method to the one outlined for the Alfven resonance is followed. 
Assuming that the inhomogeneous transition region over which density varies is sufficiently 
small, compared to the tube radius, jump conditions can be used to obtain analytical expres- 
sions for the dispersion relation. The jump conditions for the Alfven resonance are given in 
Equation d 1 6b - The corresponding jump conditions for the slow resonance at the point r = r$ 
are 



k 2 ' 9 J 



[Ir] = ~ m ~^a\~2~^TZ ( ~Z~T~2 I PT> [pt\=0, at r = r s , (19) 



where co c and <9 r po|. s are the cusp frequency and the modulus of the radial derivative of 
the density profile, evaluated at the slow resonance and c 2 is the sound speed. Note that 

the factor ( c 2^,,2 J K constant everywhere in the equilibrium. These jump conditions allow 

to obtain a more general dispersion relation that now contains the contributions from both 
Alfven and slow resonances as 



m 



D((0,k z ) = -l7l y 2 --^ : «E-% — T \ — -. (20) 

wj d r po A co 2 \c 2 + v 2 J d r p or 



To obtain an analytic expression for the damping rate of the kink mode the long-wavelength 
limit has to be considered. In terms of the physically relevant quantities, the damping time 
over the period can now be cast as 



U „ 1 fpi + P 



P (/ J a) \pi-p, 



+ (k z a) 2 ( c 2 s V 1 



cos a A m \ c 2 + v 2 . J cos a s 



(21) 



Here F is the same numerical factor as in Equation i ll St . <Xa = k{ta — a) /I, and as = 
n(rs — a) /I. The term with k~ corresponds to the contribution of the slow resonance. If 
this term is dropped and m = 1 and cos = 1 are taken Equation J2U becomes equivalent 
to Equation i ll 8t which only takes into account the Alfven resonance. 

Equation I l21b can now directly be appli ed to measure th e relative contribution of each 
resonance to the total damping. To do that lSoler et alJ |2009f ) assume r& — r$ ~ a, for sim- 
plicity, so cos ~ cos as — l - The ratio of the two terms in Equation d21b is then 

Ma _ (M 2 ( c 2 , (22) 



(x d )s m 2 \c 2 +v] 

A simple calculation shows that for the wavelength of observed filament thread oscilla- 
tions, typically in between 10~ 3 < k z a < 10~', the slow resonance is irrelevant to account 
for the kink mode damping. For instance, for m = 1 and k z a = 10~ 2 Equation d22b gives 
(Td)A/{?d)s — 10 -7 . Even with the largest ratio that can be obtained, assuming a very large 
plasma-/?, i.e., c 2 — > °°, so the factor ^ c 2^! v 2 J ~ ^ 1> gives aratio of Alfven to slow continuum 
damping times of (T d ) A / (z d ) s < 1. 
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Fig. 11 Ratio of the damping time to the period, 
Tj/P, as a function of the dimensionless wavenum- 
ber, k z a, corresponding to the kink mode for lja = 
0.2. The solid line is the full numerical solution. 
The symbols and the dashed line are the results of 
the thin boundary approximation for the Alfven and 
slow resonances. The shaded region represents the 
range of typically observed values for the wave- 
lengths in promin ence oscillations. Adapted from 
ISoleretalJ 120093) . 
rcr 5 io~ 4 io~ s w 2 w~' 10° 10' io' 

This analytical result is further confirmed bv lSoleretal.l J2009f) by performing numer- 
ical computations outside the thin tube and thin boundary approximations, by solving the 
full resistive eigenvalue problem. Figure[TT|displays the individual contribution of each res- 
onance. The slow resonance is seen to be much less efficient than the Alfven resonance. 
For the wavenumbers relevant to observed prominence oscillations, the value of Td/P due 
to the slow resonance is in between 4 and 8 orders of magnitude larger than the same ra- 
tio obtained for the Alfven resonance. On the other hand, the complete numerical solution 
(solid line) is close to the result for the Alfven resonance. As seen in Sect. I4.ll one can 
obtain values of Td/P — 3 in the relevant range of wavenumber. Note also that for short 
wavelengths (k z a > 10°) the value of Td/P increases and the efficiency of the Alfven reso- 
nance as a damping mechanism decreases. For the larger wavelengths considered to produce 
Figure [TT] k z a ~ 10 2 , both the slow and the Alfven resona nces produce simil ar (and inef- 
ficient) damping rates. The overall conclusion obtained bv lSoleretal.1 12009 j) is therefore 
that, although the plasma- jS in solar prominences is definitely non-zero, the slow resonance 
is very inefficient in damping the kink mode for typical prominence conditions and in the 
observed wavelength range. The damping times obtained with this mechanism are compara- 
ble to those due to the thermal effects discussed previously. Therefore, the resonant damping 
of transverse thread oscillations is governed by the Alfven resonance. 

5 Resonant absorption in a partially ionized filament thread 

From the results described so far resonant absorption in the Alfven continuum seems to be 
the most efficient damping mechanism for the kink mode and the only one that can produce 
the observed damping time-scales. On the other hand, the effects of partial ionization could 
also be relevant, at least for short wavelengths. The qu estion arise s on wh ether partial ion- 
ization affects the mechanism of resonant absorption. ISoler et al] l2009dl ) have integrated 
both mechanisms in a non-uniform cylindrical filament thread model in order to assess both 
analytically and numerically the combined effects of partial ionization and Alfvenic reso- 
nant absorption on the kink mode damping. Apart from the inherent relevance of this work 
in connection to the damping of prominence oscillations, this study constitutes the first of 
the kind that considers the resonant absor ption phenomenon in a partially ionized plasma. 

The filament thread model used by [Soleretal. (2009d) is an infinite straight cylinder 
with prominence-like conditions embedded in an unbounded coronal medium. Resonant 
damping is included in the model by connecting the prominence and coronal densities by 
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means of a transitional layer, with a characteristic inhomogeneity length-scale I (see expres- 
sions II II and 1121 ). The plasma properties are now also characterized by the ionization 
fraction, p.Q. The coronal plasma is assumed to be fully ionized, but the prominence material 
is only partially ionized. As with the transverse density variation, the radial behavior of the 
ionization fraction in filament threads is unknown, but a one-dimensional transverse profile, 
similar to the one used to model the equilibrium density, can be assumed. The following 
profile for the ionization fraction, p.(r), is adopted 

!p,f, if r<a — l/2 7 
fL tI {r),ifa-l/2<r<a + l/2, (23) 
p. c , if r>a + l/2, 

with 




where the filament ionization fraction, p,f, is considered a free parameter and the corona 
is assumed to be fully ionized, so p. c = 0.5. The non-uniform transitional layer of length I 
therefore connects plasma with densities in the range pf and p c and ionization degrees in 
the range jlj and p. c . As before, the one-fluid approximation is used and for simplicity the 
/3 = limit, which excludes slow waves, is considered. The quantities tj, T]c, and T)h are 
now functions of the radial direction. 



5.1 Analytical results 

In order to obtain some analytical approximations, first, only Cowling's diffusion is con- 
sidered. This basically introduces perpendicular magnetic diffusion and allows to derive an 
analytical dispersion relation for transverse oscillations, since the induction equation can be 
written in a compact form as follows 



dB, r, 



,-Vx vixBq , (25) 
dt v z A 

where = v 2 A — i(or\c is the modified Alfven speed squared dForteza et al] |2008h . The 
dispersion relation for trapped waves is obtained by imposing that the solutions are regu- 
lar at r = 0, the continuity of the radial displacement, | r , and the total pressure pertur- 
bation, pr, at the tube b oundary, and that perturbations must vanish at infinity (see, e.g., 
Ed win and Roberts| [l983). The dispersion relation reads D(co,k z ) = 0, with 

D{coM = j 2 n \*' mi r\ - ? mf 7 /:f fa l (26) 
p c (0)2 - £?r A 2 c ) K m (n c a) pf f m 2 _ k 2p2\ J m (m f a) 

where /,„ and K m are the Bessel function and the modified Bessel function of the first kind, 
respectively, and the quantities m / and n c are given by 



*/ = ^-f2 "c = ^jr 2 '- (27) 

A/ ^Ac 

As shown in Section |4] an analytical dispersion relation in the case I /a 7^ can be 
obtained in the thin boundary approximation, by making use of the jump conditions for the 
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radial displacement and the total pressure perturbation at the Alfven resonance (expressions 
1161 ). This dispersion relation is 

m 2 Ir 2 

D(a,k z ) = -in A , . (28) 

<WP0|r A 

Note that Equation ( 128b is formally identical to Equation J 17b . but now D(co,k z ) is de- 
fined in Equation J26b . By combinin g the thin tube and thin boundary approximations, and 
neglecting terms proportional to fc?, Soler et al] |2009d) arrive at the following short-hand 



expression for the damping time over the period 



P % 



m | L\ ( Pf-PA + 2{pfflcf + Pcflcc)k z a 



aJ\p f + Pc J / 2p/ (p / + p , 



(29) 



with f\cc.f = 1c/vAc,/a the filament and coronal Cowling's diffusivities in dimensionless 
form. We see that the term related to the resonant damping is independent of the value of 
Cowling's diffusivity and, therefore, of the ioni zation degree, and take s the same form as in 
a fully ionized plasma, s ee e.g . Equation (77) in lGoossens et al.l dl992h and Equation (56) in 
iRuderman and Roberts! J2002h . On the other hand, the second term related to the damping 
by Cowling's diffusion is proportional to k z , so its influence in the long-wavelength limit is 
expected to be small. To perform a simple calculation, consider the case m = 1, k z a = 10 -2 , 
and I /a = 0.2. This results in Ta/P ~ 3.18 for a fully ionized thread (fif = 0.5), and T^/P w 
3.16 for an almost neutral thread (p.f = 0.95). The ratio T^/P depends very slightly on the 
ionization degree, suggesting that resonant absorption dominates over Cowling's diffusion. 

The relative importance of the two mechanisms can be estimated by taking the ratio of 
the two terms on the right-hand side of Equation i 



(^)ra _ / 2 (Pf + Pc) ( Pffkf + Pclcc \ k z a 
(T d ) c V Pf V Pf-Pc ) m{l/aY 

with (Td) RA and (Td) c the damping times due to resonant absorption and Cowling's dif- 
fusion, respectively. This last expression can be further simplified by considering that in 
filament threads py 2> p c and f/c/ 2> r/c c , so that 

k, V2t]c/ , , ■ (31) 
(Td)c m(l/a) 

The efficiency of Cowling's diffusion with respect to that of resonant absorption increases 
with k z a and p.f (through f/c/)- Considering the same parameters as before, one obtains 
(T d ) RA / (T d ) c w 2 x 10 -8 for fif = 0.5, and (Td)^ / (T d ) c « 6 x 10~ 3 for fLf = 0.95, mean- 
ing that resonant absorption is much more efficient than Cowling's diffusion. From Equa- 
tion OH it is also possible to estimate the wavenumber for which Cowling's diffusion be- 
comes more important than resonant absorption, by setting (Td) RA / (^d)c ~ 1= as 

m(l/a) 

kMK X ■ (32) 
V2f] C f 

Considering again the same parameters, Equation 432b gives k,a w 5 x 10 5 for fi f = 0.5, 
and k z a w 1.7 for p.f = 0.95. Note however that Equation ( 13 lb is only valid for k z a <?C 1, 
so resonant absorption is expected to be the dominant damping mechanism in the long- 
wavelength regime, even for an almost neutral filament plasma. 
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Fig. 12 Ratio of the damping time to the period of the kink mode as a function of k z a corresponding to 
a thread without transitional layer, i.e., 1/a = 0. (a) Results for a = 100 km considering different ionization 
degrees: p.f — 0.5 (dotted line), flf = 0.6 (dashed line), p.f = 0.8 (solid line), and fif = 0.95 (dash-dotted line). 
Symbols are the approximate solution given by solving Equation (26) for flf = 0.8. (b) Results for flf = 0.8 
considering different thread widths: a = 100 km (solid line), a = 50 km (dotted line), and a = 200 km (dashed 
line). The sha ded zone correspond s to the range of typically observed wavelengths of prominence oscillations. 
Adapted from lSoler eta?] <2009dh . 



5.2 Numerical results 

The analytical estimates described above are verified and extended bv lSoler et alJ d2009dh by 
numerically solving the full eigenvalue problem. Computations include now Hall's diffusion 
in addition to ohmic and Cowling's dissipation. 

First, a homogeneous filament thread without transitional layer (I /a = 0) is considered. 
This is the same case presented in Section 13.21 with the addition of Hall's term in the in- 
duction equation. Figure [T2"a displays the obtained damping rate for different ionization 
degrees. In agreement with the results displayed for the kink mode in Figure [8] i^jP has a 
maximum which corresponds to the transition between the ohmic-dominated regime, which 
is almost independent of the ionization degree, to the region where Cowling's diffusion is 
more relevant and the ionization degree has a significant influence. The approximate solu- 
tion obtained by solving Equation d26b for a given value of p. agrees very well with the 
numerical solution in the region where Cowling's diffusion dominates, while it significantly 
diverges from the numerical solution in the region where ohmic diffusion is relevant. Within 
the range of typically reported wavelengths, is between 1 and 2 orders of magnitude 
larger than the measured values, so neither ohmic nor Cowling's diffusion can account for 
the observed damping time. Figure IT2b shows that the smaller the thread radius, the more 
rapidly the kink wave is attenuated. In addition, the critical wavenumbers are shifted when 
varying the thickness of the threads and the wavenumber range for which the kink wave 
propagates is wider for thick threads than for thin threads. The critical wavenumbers are far 
from the relevant values of k z a for the observed thread widths, and so they should not affect 
the kink wave propagation in realistic filament threads. 

Next, the inhomogeneous case (/ /a ^ 0) is analyzed, with the inclusion of resonant wave 
coupling in the Alfven continuum. Figure [T3k displays some relevant differences with re- 
spect to the homogeneous case (I /a = 0). First, the damping time is dramatically reduced for 
intermediate values of k z a including the region of typically observed wavelengths. In this re- 
gion, the ratio T&/P becomes smaller as I /a is increased, a behavior consistent with damping 
by resonant absorption. The ratio i&jP is independent of the thickness of the layer, and for 
large k z a coincides with the solution in the homogeneous case. The cause of this behavior 
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Fig. 13 Ratio of the damping time to the period of the kink mode as a function of k z a corresponding to a 
thread with an inhomogeneous transitional layer, (a) Results for (if = 0.8 considering different transitional 
layer widths: l/a = (dotted line), l/a = 0.1 (dashed line), l/a = 0.2 (solid line), and l/a = 0.4 (dash-dotted 
line). Symbols are the solution in the TB approximation given by solving Equation j28t for l/a = 0.2. (b) 
Results for l/a = 0.2 considering different ionization degrees: (if = 0.5 (dotted line), (if = 0.6 (dashed line), 
(if = 0.8 (sol id line), and (if = 0.95 (dash-dotted line). In both panels we have considered a = 100 km. 
Adapted from lSoler etafl <2009dh . 




Fig. 14 Ratio of the damping time to the period of 
the kink mode as a function of k z a corresponding to a 
thread with a = 100 km and l/a = 0.2. The different 
linestyles represent the results for: partially ionized 
thread with fit = 0.8 considering all the terms in 
the induction equation (solid line), partially ionized 
thread with (if = 0.8 neglecting Hall's term (sym- 
bols) , and fully ionized th read (dotted line). Adapted 
from lSoler eta?] <2009d) . 



is that perturbations are essentially confined within the homogeneous part of the thread for 
large k z a, and therefore the kink mode is mainly governed by the internal plasma conditions. 
On the other hand, the solution for small k z a is completely different when l/a =^0. The in- 
clusion of the inhomogeneous transitional layer removes the smaller critical wavenumber 
and consequently the kink mode exists for very small values of k z a. Figure[T3h also shows a 
very good agreement between the numerical and the analytical solutions (obtained by solv- 
ing the dispersion relation 1281 ). for wavenumbers above k z a ~ 10~ 4 , and a poor agreement 
in the range for which ohmic diffusion dominates, below k z a ~ 10~ 4 . To understand this one 
has to bear in mind that dispersion relation d28b takes into account the effects of resonant 
absorption and Cowling's diffusion, but not the influence of ohmic diffusion. As Figure IT3b 
shows, the ionization degree is only relevant for large wavenumbers, where the damping rate 
significantly depends on the ionization fraction. 

Finally, Figure[l4]displays the ranges of k z a where Cowling's and Hall's diffusion dom- 
inate. As expected, Hall's diffusion is irrelevant in the whole studied range of k z a, while 
Cowling's diffusion dominates the damping for large k z a. In the whole range of relevant 
wavelengths, resonant absorption is the most efficient damping mechanism, and the damping 
time is independent of the ionization degree as predicted by the analytical result (Eq. 1291 ). 
On the contrary, ohmic diffusion dominates for very small k z a. In that region, the damping 
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time related to Ohm's dissipation becomes even smaller than that due to resonant absorption, 
meaning that the kink wave is mainly damped by ohmic diffusion. 



6 Application to prominence seismology 

Solar atmospheric seismology aims to determine physical parameters in magnetic and plasma 
structures that are difficult to measure by direct means. It is a remote diagnostics method 
that combines observations of oscillations and waves in magnetic structures, together with 
theoretical results from the analysis of oscillatory properties of models of those structures. 
The philosophy behind this discipline is akin to that of Earth seismology, the sounding 
of the Earth interior using seismic waves, and helio-se i smolo gy, t he acoustic diagnosti c 



01 the fc,arth interior using seismic waves, and helio-seismology, the acoustic diagnostic 
of the solar interior. It was first sugge sted by lUchidal dl970h and iRoberts et all dl984) 



in the coronal context, and by iTandberg-Hanssenl ( 1 19951) in a prominence context. The 
last years increase in the number and quality of high resolution observations has lead to 
its rapid development. In the context of coronal loop oscillations, recent applications of 
this technique have allowed t he estimation and/or re striction of parameters such as the 
magnetic field strength dNakariakov and Ofmanll200l|). the Alfven speed in coronal loops 
dZaqar ashvili 2003; Arreg uTet alJ2007l : lGoossens et al.ll2008h . t he transversal density struc- 
turing Jverwichte et alJ2006h . or the coronal density scale height dAndries et alj2005llVerth et al.1 
2008). Its application to prominences is less developed. Some recent results of the MHD 



prominence seismology technique are shown. 



6.1 Seismology using the period of filament thread oscillations 

The first prominence seismology application u sing Hinode observat ions of flowing and 
transversely oscillating threads was prese nted bylj erradas et al. (2008), using observations 
obtained in an active region filament by lOkamoto et alj d2007l) . The observations show a 
number of threads that flow follo wing a path pa r allel t o the photosphere while they are os- 
cillating in the vertical direction. Terra das et al.l d2008l) interpret these oscillations in terms 
of the kink mode of a magnetic flux tube. By using previous theoretical res ults from a norma l 
mod e analysis in a two-dimensiona l piecewise filament thread model by iDfaz et al.l |2002) 
and iDvmova and Rudermanl d2005l) , these authors find that, although it is not possible to 
univocally determine the physical parameters of interest, a one-to-one relation between the 
thread Alfven speed and the coronal Alfven speed can be established. This relation comes 
in the form of a number of curves relating the two Alfven speeds for different values of the 
length of the magnetic flux tube and the density contrast between the filament and coronal 
plasma. The obtained curves have an asymptotic behavior for large values of the density con- 
trast, typical of filament to coronal plasmas, and hence a lower limit for the thread Alfven 
spe ed can be obt ained. Further details on this study can be found in iTerradas et all (2008) 
and lOlivedd2009l) . 

A recent application of the prominence seismology techni que, using the p eriod of ob- 
served filament thread transverse oscillations can be found in lLin et alld2009l) . These au- 
thors find observational evidence about swaying motions of individual filament threads from 
high resolution observations obtained with the Swedish 1-m Solar Telescope in La Palma. 
The presence o f waves propagating along individual threads was already evident in, e.g., 
iLin et alj d2007l) . However, the fact that line-of-sight oscillations are observed in promi- 
nences beyond the limb, as well as in filaments against the disk, suggests that the planes 
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Fig. 15 (a): Ratio c?/vjL (solid line) as a function of the density contrast, c. The dotted line corresponds to 
the value of the ratio cf/v^ for c — > °°. (b): Magnetic field strength as a function of the internal density, p,, 
corresponding to some selected threads: thread 1 (solid line), thread 3 (dotted line), thread 5 (dashed line), 
and thread 7 (dash-dotted line). Adapted from Lin et al. 1 200^). 



of the oscillation ma y acquire various orientations relative to the local solar reference sys- 
tem. For this reason, iLin et al.l ( l2009h combine simultaneous recordings of motions in the 
line-of-sight and in the plane of the sky, which leads to information about the orientation of 
the oscillatory plane in each case. Periodic oscillatory signals are obtained in a number of 
threads, that are next fitted to sine curves, from which the period and the amplitude of the 
waves are derived. The presence of different cuts along the structures allow llin et alj 12009) 
to obtain the phase difference between the fitted curves, which can be used to measure the 
phase velocities of the waves. The overall periods and mean velocity amplitudes that are ob- 
tained correspond to short period, P ~ 3.6 s, and small amplitude ~ 2 km s oscillations. 
The information obtained from these H a filtergrams in the plane of the sky is combined with 
H a Dopplergrams, which allow to detect oscillations in the line-of-sight direction. By com- 
bining the observed oscillations in the two orthogonal directions the full vectors are derived, 
which show that the oscillatory planes are close to the vertical. 

ILin et al l d2009h interpret the observed swaying thread oscillations as kink MHD waves 
supported by the thread body. By assuming the classic one-dimensional, straight, flux tube 
model a comparison between the observed wave properties and the theoretical prediction is 
performed on order to obtain the physical parameters of interest, namely the Alfven speed 
and the magnetic field strength in the studied threads. To this end c k = (o/k z with CO defined 
by Equation ( 1 14b is used as an approximation to the kink speed, which can directly be 
associated to the measured phase velocity of the observed disturbances. Figure [T5k shows 
that, for the large density contrasts expected for filament threads, the curve that relates the 
kink speed to the Alfven speed becomes flat, and so the ratio is almost independent of 
the density contrast. This allows to simplify the kink speed to cj ~ V2va, and therefore 
obtain the thread Alfven speed thro ugh vai — Khase / "\/2. The obtained values for a set of 10 
threads can be found in Table 2 in llin et all (2009). Once the Alfven speed in each thread 
is determined, the magnetic field strength can be computed, if a given internal density is 
assumed. For a typical value of p, = 5 x 10~ n kg m~ 3 magnetic field strengths in between 
0.9 -3.5 G are obtained, for the analyzed events. If the unknown density is allowed to vary 
in a range of plausible values, the derived magnetic field strength also changes accordingly, 
as can be seen i n Figure [15b. The important conclusion that we extract from the analysis by 
ILin et alj (2009) is that prominence seismology is possible and works well, provided high 
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Fig. 16 {Left): Analytic inversion of physical parameters in the (c, 1/a, Va/) space for a filament thread os- 
cillation with P = 3 min, = 9 min, and a wavelength of A = 3000 km (see e.g. Lin et al. 2007). (Right): 
Magnetic field strength as a function of the density contrast and transverse inhomogeneity length-scale, de- 
rived from the analytic inversion for a coronal density of p c = 2.5 x 10~ 13 kg m~ 3 . 

resolution observations are available. The derived plasma parameters show rather different 
values for different threads, which indicates that the plasma parameters may be varying 
along individual threads belonging to the same filament. This however does not come as a 
surprise in view of the highly inhomogeneous nature of these objects. 

6.2 Seismology using the period and damping of filament thread oscillations 

A feature clearly observed by llin etal ] d2009h is that the amplitudes of the waves pass- 
ing through two different cuts along a thread are notably different. Apparent changes can 
be due to damping of the waves in addition to the noise in the data. Among the different 
damping mechanisms that are described in this paper, resonant absorption in the Alfven 
continuum seems a very plausible mechanism and is open to direct application for filament 
thread seismology, using the dampi ng as an additio n al source of informatio n. This additional 
information has been used bv lArregui et al.l d2007t) : lGoossens et al. N2008h in the context of 
transversely oscillating coronal loops and its application to filament threads is straightfor- 
ward. 

Th e analytical and numerical inversion schemes bv lArregui et al.ll2007l) and lGoossens et all 
(2008) make use of the simple idea that it is the same magnetic structure, whose equilibrium 
conditions we are interested to assess, that is oscillating with a given period and undergoing 
a given damping rate. By computing the kink normal mode frequency and damping time as 
a function of the relevant equilibrium parameters for a one-dimensional model, the period 
and damping rate have the following dependencies 



where Ta; the internal Alfven travel time. In the case of coronal loop oscillations, an estimate 
for k z can be obtained directly from the length of the loop and the fact that the fundamental 



P = P(k z ,c,l/a,X Ai ), 




(33) 
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kink mode wavelength is twice this quantity. For filament threads, the wavelength of oscil- 
lations needs to be measured. Relations d33b indicate that, if no assumption is made in any 
of the physical parameters of interest, we have two observed quantities, period and damp- 
ing time, and three unknowns, density contrast, transverse inhomogeneity length-scale, and 
Alfven travel time (or conversely Alfven speed through the relation = L/tai). There 
are therefore infinite different equilibrium models that can equally well explain the obser- 
vations. These valid equilibrium models are displayed in Figure [T6k. wher e the analytical 
algebr aic expressions in the thin tube and thin boundary approximations by Goo ssens et al.l 
(2008) have been used to invert the problem. It can be appreciated that, even if an infi- 
nite number of solutions is obtained, they define a rather constrained range of values for 
the thread Alfven speed. Also, and because of the insensitiveness of the damping rate with 
density contrast, for the typically large values of this parameter in prominence plasmas, the 
obtained solution curve displays an asymptotic behavior for large values of c. This allows 
us to obtain precise estimates for the filament thread Alfven speed, Vai — 12 km s , and the 
transverse inhomogeneity length scale, I /a ~ 0.16. Note that these asymptotic values can 
directly be obtained by inverting Equations dl5b and d 1 8b for the period and the damping 
rate in the limit c — > °°. The computation of the magnetic field strength from the obtained 
seismological curve requires the assumption of a particular value for either the filament or 
the coronal density. The resulting curve for a typical coronal density is shown in Figure[T6b. 
As can be seen, precise values of the magnetic field strength cannot be obtained, unless the 
density contrast is accurately known. 



7 Open Issues 

Solar prominences are, probably, the most complex structures in the solar corona. A full un- 
derstanding of their formation, magnetic structure and disappearance has not been reached 
yet, and a lot of physical effects remain to be included in prominence models. For this reason, 
theoretical models set up to interpret small-amplitude oscillations and their damping are still 
poor. High-resolution observations of filaments suggest that they are made of threads whose 
thickness is at the limit of the available spatial resolution. Then, one may wonder whether 
future improvements of the spatial resolution will provide with thinner and thinner threads 
or, on the contrary, there is a lower limit for thickness and we will able to determine it in 
the near future. The presence of these long and thin threads together with the place where 
they are anchored and the presence of flows along them suggest that they are thin flux tubes 
filled with continuous or discontinuous cool material. This cool material is probably sub- 
ject to cooling, heating, ionization, recombination, motions, etc. which, altogether, makes 
a theoretical treatment very difficult. For instance, in the case of the considered thermal 
mechanisms, up to now only optically thin radiation has been considered while, probably, 
the consideration of optically thick effects would be more realistic; the prominence heating 
mechanisms usually taken into account are tentative and "ad hoc" while true prominence 
heating processes are still deeply unknown. An important step ahead would be to couple 
radiative transfer with magnetohydrodynamic waves as a means to establish a relationship 
between velocity, density, magnetic field and temperature perturbations, and the observed 
signatures of oscillations like spectral line shift, width and intensity. Partial ionization is an- 
other topic of interest for prominence oscillations since, apart from influencing the behavior 
of magnetohydrodynamic waves, it poses an important problem for prominence equilib- 
rium models when gravity is taken into account and which is: How non-ionized prominence 
material is kept in the magnetic field?. Another issue which still remains a mystery is the 
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triggering mechanism of small-amplitude oscillations. In the case of large-amplitude oscil- 
lations, observations provide with information about the exciting mechanism, however, the 
available observations of small-amplitude oscillations do not show any signature of the ex- 
citing mechanism. Are these oscillations of chromospheric or photospheric origin? Are they 
generated inside prominence magnetic structure by small reconnection events? Are they pro- 
duced by weak external disturbances coming from far away in the solar atmosphere? The 
changing physical conditions of prominence plasmas suggest that for an in-depth theoretical 
study of prominence oscillations more complex models together with numerical simula- 
tions are needed. Therefore, and as a step ahead, in the next future numerical studies of the 
time evolution of magnetohydrodynamic waves in partially ionized flowing inhomogeneous 
prominence plasmas, subject to different physical processes such as ionization, recombina- 
tion, etc., should be undertaken. However, a full three-dimensional dynamical prominence 
model involving magnetic equilibrium, radiative transfer, etc. whose oscillatory behavior 
could be studied seems to be still far away in the future. 
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